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Abstract 

Background: The microbiota of the mammalian gastrointestinal (Gl) tract consists of diverse populations of 
commensal bacteria that interact with host physiological function. Dysregulating these populations, through 
exogenous means such as antibiotics or dietary changes, can have adverse consequences on the health of the 
host. Studies from laboratories such as ours have demonstrated that exposure to psychological stressors disrupts 
the population profile of intestinal microbiota. To date, such studies have primarily focused on prolonged stressors 
(repeated across several days) and have assessed fecal bacterial populations. It is not known whether shorter 
stressors can also impact the microbiota, and whether colonic mucosa-associated populations can also be affected. 
The mucosa-associated microbiota exist in close proximity to elements of the host immune system and the two are 
tightly interrelated. Therefore, alterations in these populations should be emphasized. Additionally, stressors can 
induce differential responses in anxiety-like behavior and corticosterone outputs in variant strains of mice. Thus, 
whether stressor exposure can have contrasting effects on the colonic microbiota in inbred C57BL/6 mice and 
outbred CD-I mice was also examined. 

Results: In the present study, we used high throughput pyrosequencing to assess the effects of a single 2-hour 
exposure to a social stressor, called social disruption (SDR), on colonic mucosa-associated microbial profiles of 
C57BL/6 mice. The data indicate that exposure to the stressor significantly changed the community profile and 
significantly reduced the relative proportions of two genera and one family of highly abundant intestinal bacteria, 
including the genus Lactobacillus. This finding was confirmed using a quantitative real-time polymerase chain 
reaction (qPCR) technique. The use of qPCR also identified mouse strain-specific differences in bacterial abundances. 
L reuteri, an immunomodulatory species, was decreased in stressor-exposed CD-I mice, but not C57BL/6 mice. 

Conclusions: These data illustrate that stressor exposure can affect microbial populations, including the lactobacilli, 
that are closely associated with the colonic mucosa. Because the lactobacilli can have beneficial effects on human 
health, stressor-induced reductions of their population could have important health implications. 
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Background 

The human body is colonized by an enormous array of 
microbes collectively referred to as the microbiota. It is es- 
timated that there are approximately ten times more bac- 
terial cells than there are human cells in the human body 
with microbiota being found in various receptive niches 
such as on the skin, in the oral and respiratory tracts, re- 
productive tract, and most numerously, in the gastrointes- 
tinal (GI) tract [1,2]. In the GI tract, the microbiota form a 
community that fully interacts with one another and with 
the host. The GI microbiota have multiple beneficial ef- 
fects on their hosts [3-5] and the structure of the GI 
community can impact these functions. Of the benefi- 
cial interactions between microbiota and host, those 
that involve the immune system are some of the most 
intensively studied, and several lines of experiments 
have demonstrated that the GI microbiota and the host 
immune system are deeply intertwined in each other s 
development [6]. For example, germ-free mice, which 
have never been colonized by any type of microbe, have 
diminished immune responses compared to mice colo- 
nized with a healthy microbiota [7,8]. Mice deficient in 
Nod2, a gene that encodes for a receptor that is involved 
in immune recognition of bacterial muramyl dipeptide, 
develop a unique microbiota that can lead to colonic in- 
flammation when transmitted to wild-type mice [9,10]. 
It is now well recognized that the GI microbiota have 
strong effects on immunoregulation, and immune system 
activity, in turn, helps to shape the GI microbiota [6-9]. 

It is estimated that over 500 different species of bacteria 
can colonize the GI tract [11]. Despite this enormous 
species-level diversity, only a small number of bacterial 
phyla are represented. The vast majority of bacteria in 
the human GI tracts belong to the phyla Firmicutes 
and BacteroideteSy with bacteria in the phyla Actinobac- 
teria and Proteobacteria also comprising a smaller portion 
of the overall microbiota as well [12]. The GI microbiota 
of the murine GI tract is similar to that of humans, with 
the majority of microbes belonging to the major phyla, 
Firmicutes and Bacteroidetes [13-15]. These populations 
shift slowly over time, but their general stability is import- 
ant to the health of the host. Abrupt changes to the gut 
microbiota have been shown to potentially lead to serious 
negative host health outcomes, including diarrhea, oppor- 
tunistic infections, and obesity [16-18]. Changes to the 
microbiota can be caused by factors such as antibiotic use 
and severe enteric infection. Data from this laboratory, as 
well as others, indicate that exposure to either physical or 
psychological stressors can also alter intestinal microbe 
profiles. Stressor exposure early in life has been demon- 
strated to alter the types and abundance of bacteria found 
in the intestines. Separating infant monkeys from their 
mothers to induce a physiological stress response resulted 
in a significant reduction in the number of total lactobacilli 



that could be cultured from the stool [19]. Reductions in 
lactobacilli are meaningful as certain species, including L, 
reuteriy are involved in immunomodulation [20,21]. In 
rats, separating the pups from their mothers during the 
first 14 days of life led to an altered GI microbiome [22]. 
Stressor exposure during adulthood can also impact the 
stability of the intestinal microbiota. For example, expos- 
ing adult mice to a prolonged restraint stressor was shown 
to significantly alter microbial profiles in the cecal con- 
tents [23]. 

Similar findings have been associated with a social stres- 
sor called social disruption (SDR) [15] that involves re- 
peated social defeat as a result of inter-male aggression 
over a 2 hr period. When repeated over 6 consecutive 
nights, this stressor induces a physiological stress response 
marked by the activation of the hypothalamic pituitary 
adrenal (HPA) axis and the sympathetic nervous system 
[24]. The study by Bailey et al. demonstrated that this 
week-long exposure to the stressor can alter the cecal 
luminal microbiome in outbred CD-I mice [15]. How- 
ever, it is not known whether a single 2 hr exposure to 
the stressor could induce similar alterations, or whether 
other mouse strains, such as inbred C57BL/6 mice that 
are widely used in infectious disease research, are affected 
by stressor exposure. Knowing whether short-lasting stressors 
can impact the microbiota has translational importance 
particularly for the inflammatory bowel diseases (IBD) 
and irritable bowel syndrome (IBS). Patients with IBD 
often report stressful periods preceding symptom flares 
[25]. It is not clear whether acute stress can exacerbate 
IBD symptomatology, but it is recognized that acute 
stressor exposure is associated with increased visceral 
sensitivity in IBS [26] . Because alterations in gut micro- 
biota community structure are evident in IBD and IBS 
[27,28], and because these alterations are thought to 
possibly contribute to these diseases [29,30], we deter- 
mined whether a short-lasting stressor was sufficient to 
impact gut microbiota community structure. 

Previous studies assessing the effects of stress on the 
microbiota have relied on assessment of intestinal lumen 
or fecal populations [15,23,31]. Studies in healthy indi- 
viduals as weU as patients with IBS, IBD, or hepatic en- 
cephalopathy demonstrate that luminal/fecal microbiota 
are significantly different than mucosa-associated micro- 
biota [32,33]. Because shifts in microbial populations that 
are in close proximity to the intestinal epithelium are 
thought to have the strongest effects on host immunity 
[34], an additional aspect of this study was to determine 
whether stressor exposure could impact tissue-associated 
microbial community profiles. 

Stressor-induced changes in lactobacilli have reliably been 
found in laboratory animals as well as in humans [19,35,36]. 
However, stressor-induced changes in other taxa have 
not been widely reported. Thus, 454 pyrosequencing 
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was performed on the colonic tissue of inbred C57BL/6 
mice to determine whether a single 2 hr stressor exposure 
could change colonic tissue-associated microbial commu- 
nity profiles. Quantitative PGR (qPCR) was then used to 
determine whether stressor-induced alterations in relative 
abundances evident with the 454 pyrosequencing were 
also evident as alterations in absolute abundance in both 
inbred and outbred mice. 



Results 

Social stress affects the community structure of the 
colonic mucosa-associated microbiota 

We analyzed the effect of a short term SDR stressor on 
the colonic mucosa-associated microbiota using 16S 
rRNA gene pyrosequencing of samples from both SDR 
C57BL/6 stressor-exposed mice and non-stressed HCC 
control mice. After normalization of the datasets by rar- 
efaction, we observed no statistical difference in the total 
number of operational taxonomic units (OTU) between 
the SDR and control mice (data not shown). Additional 
analysis of the Shannon, equitability and Chaol alpha di- 
versity metrics showed no statistical difference between the 
two groups (Figure lA-C). Thus, a single two-hour cycle 
of exposure to the social stressor did not affect the alpha 
diversity of the mucosa-associated microbiota. 

While we did not observe any statistically significant 
differences in the alpha diversity of the microbiota be- 
tween the SDR and control group, beta diversity analyses 
did reveal differences in composition. A principal coordi- 
nates analysis (PCoA) of the UniFrac distances between 
samples showed that they clustered into separate groups 
according to treatment (Figure 2A). An unweighted pair 
group with arithmetic mean (UPGMA) hierarchical clus- 
tering dendrogram confirmed this clustering (Figure 2B). 
Statistical analysis of the UniFrac distances using the ana- 
lysis of similarity statistic (ANOSIM) indicated that colonic 
mucosal microbiota from SDR stressor-exposed mice were 



significantly different from those of the home cage control 
(HCC) mice (p < 0.05). 

Populations of the genus Lactobacillus and L reuteri are 
reduced in mice exposed to a social stressor 

The taxonomic profile of the microbiota at the phylum 
level showed no significant differences between the two 
treatment groups (Figure 3). Analysis of relative abun- 
dances at lower taxonomic levels showed a reduction in 
the family Porphyromonadaceae in SDR stressor-exposed 
mice compared to non-stressed HCC mice (p < 0.01) 
(Table 1). Exposure to the SDR stressor also reduced the 
relative abundance of bacteria in the family Lactobacilla- 
ceae (p < .05) (Table 1). This reduction in the Lactobacilla- 
ceae was reflected by a reduction in the relative abundance 
of bacteria in the genus Lactobacillus (p < .05) (Table 2). 
In addition to the reduction in the abundance of lactoba- 
cilli, it was evident that exposure to the SDR stressor sig- 
nificantly reduced the relative abundance of the genus 
Parabacteroides (p < .01), as weU as an unclassified group 
within the phylum Firmicutes (p < .05) and an unclassified 
group within the class Bacilli (p < .05) (Table 2). 

Many bacteria in the genus Lactobacillus have been 
shown to have immunomodulatory ftinctions in the colon, 
and can be used as a probiotic to treat inflammation. Thus, 
we used qPCR to quantify the effect of the SDR stressor 
on the abundance of this group of bacteria [20,21]. When 
the a priori hypothesis that repeated administration of the 
SDR stressor would reduce the abundance of lactobacilli 
was tested, it was determined that the absolute abundance 
of Lactobacillus spp. was significantly lower (p < .05) after 
6 cycles of the stressor in comparison to baseline levels in 
C57BL/6 mice (Figure 4A). This reduction in lactobacilli 
abundance was also observed in outbred CD-I mice 
(p < .05), with the largest reduction in mean levels of 
lactobacilli also occurring after 6 days of the SDR stressor 
(Figure 4B). In addition, qPCR analyses also revealed that 
exposure to the SDR stressor significantly reduced the 
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Figure 1 Exposure to the SDR stressor did not impact alpha-diversity. Mice were exposed to the SDR stressor, or were left undisturbed as 
non-stressed Home Cage controls (HCC control). After 454 pyrosequencing, three measures of community alpha-diversity were calculated using 
QIIME. A. Equitability Index, B. Shannon Diversity Index, and C. Chaol Rarefaction Measurement were unaffected by exposure to the SDR stressor. 
Data are from n = 5 mice per group. 
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Principle Coordinate Analysis of Effect of Social Disruption Stress on Colonic Microbiota 
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Figure 2 Exposure to the SDR stressor significantly changes beta-diversity. After 454 pyrosequencing, UniFrac distances were calculated 
using QIIME and plotted on a 2-D principal coordinates graph (A). Exposure to the SDR stressor shifted the colonic microbiota, with 4 of 5 SDR 
stressor samples clustering apart from HCC controls. UniFrac distances were also used for construction of an unweighted pair group method by 
arithmetic mean (UPGMA) dendrogram, in which 4 of 5 SDR stressor samples clustered together (B). The clustering in Figure 2A and 2B was 
statistically significant (p < 0.05) after performing the ANOSIM test to measure differences between distance matrices of SDR stressor and HCC 
control groups. Data are from n = 5 mice per group. 
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abundance of the immunomodulatory species L reuteri 
over the course of six cycles of SDR (p < .05). This differ- 
ence was only evident in outbred CD-I mice (Figure 4C). 
L. reuteri levels were below the qPCR limit of detection of 
4.5 copies/gram of wet tissue (logio) in inbred C57BL/6 
mice. 

Relative bacterial levels differ by mouse strain 

In order to confirm that stressor exposure was reducing 
the absolute abundance of other bacterial groups that 
were reduced in relative abundance in the pyrosequenc- 
ing analysis, qPCR was performed on colonic tissues 
from both inbred C57BL/6 mice and outbred CDl 
mice. Primers targeting the 16 s ribosomal RNA genes 



of Parabacteroides distasonis, a member of the Parabac- 
teroides genus, and Bacteroides-Prevotella-Porphyromonas 
were used. Stressor exposure did not significantly affect 
the absolute abundance of either P, distasonis or Bac- 
teroides-Prevotella-Porphyromonas genera (Table 3). In 
addition, levels of Bacteroides-Prevotella-Porphyromonas 
group and P, distasonis were similar in the two strains of 
mice. 

Colonic inflammatory cytokine mRNA levels are not 
affected by social stressor exposure 

The colonic mRNA levels for IL-1|3, TNF-a, and iNOS 
were unaffected by exposure to the SDR stressor, sug- 
gesting that exposing C57BL/6 and CD-I mice to the 
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Figure 3 The relative abundances of mucosal-associated bacterial phyla were unaffected by exposure to the SDR stressor. Data are the 
mean+/-S.E. of the relative abundances calculated from 454 pyrosequencing data and are from n = 5 mice per group. 
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Table 1 Top 10 most abundant colonic tissue-associated 
bacterial families 





HCC 


SDR 


Lachnospiraceae 


38.71 ±5.22 


54.10 ±5.81 


Clostridiales; Other 


9.01 ±0.56 


10.1 4 ±0.96 


Bacteroidaceae 


10.77 ± 1.49 


6.67 ±2.11 


Loctobocilloceoe 


13.66 ±4.81 


3.13 ±100* 


Ruminococcoceoe 


5.22 ±0.36 


4.35 ± 0.62 


Acetobacteraceae 


0.001 ±0.004 


8.86 ± 8.83 


Clostridioceoe 


6.34 ±3.62 


0.92 ± 0.78 


Porphyromonodoceae 


2.67 ± 0.20 


1.59 ±0.25** 


Bacteria; Other 


2.21 ±0.25 


1.84 ±0.46 


Peptostreptococcaceae 


1.81 ±0.35 


0.94 ± 0.40 


Data are the mean ± standard error. *p < .05 vs. HCC; ''''p < .01 vs. HCC. 

SDR stressor did not result in detectable increases in co- 
lonic cytokines or inflammatory mediators (Table 4). 



Discussion 

Exposure to physical and psychological stressors has 
been shown to impact the gut microbiota of both labora- 
tory animals and humans [19,22,35]. However, the ma- 
jority of the previous studies have utilized repeated and 
prolonged stressors, and have assessed microbial popula- 
tions in the lumen of the intestines or present in the 
fecal matter. Whether stressor exposure has different ef- 
fects on the microbiota of different strains of mice is 
also poorly understood. Because gut microbes that ad- 
here to the colonic mucosa can have different effects on 
the host [34], this study assessed whether stressor expos- 
ure could alter the community structure of mucosa- 
associated microbes. This study demonstrated that as lit- 
tle as 2 hrs of stressor exposure is enough to signifi- 
cantly change the structure of the microbiota associated 
with the colonic mucosa. This effect was not manifest as 
alterations in alpha diversity, but rather as alterations in 
beta diversity. The ANOSIM distance matrix analysis 
and cluster analysis based on unweighted UniFrac demon- 
strated that microbiota within the samples from stressor- 
exposed mice were significantly different than those from 
the non-stressed HCC control mice. In addition to signifi- 
cantly changing beta diversity, exposure to the stressor 
significantly changed the relative abundance of 2 bacterial 
genera, namely Parabacteroides and Lactobacillus ^ and 
one bacterial family, Porphyromonadaceae, 

This study targeted mucosa-associated populations, 
but true stratification of luminal and mucosal popula- 
tions does not exist. There is much crossover between 
microbes often associated with the lumen and those that 
can adhere to the mucus layer of the GI tract, as the 
former becomes trapped in the mucus layer and the lat- 
ter sheds into the lumen. Steps were taken to remove 



Table 2 Top 40 most abundant tissue-associated bacterial 
genera 





HCC 




SDR 


Lachnospiraceae; Other 


37.61 : 


± 5.23 


53.02 ± 5.74 


Clostridiales; Other; Other 


9.01 ± 


0.56 


10.14 ± 0.96 


Bacteroides spp. 


10.77 : 


± 0.56 


6.67 ± 2.1 1 


Lactobacillus spp 


13.53 : 


±4.76 


3.11 ± 1.00* 


Rosemonas spp 


0.00 ± 


0.00 


8.83 ± 8.83 


Ruminococcaceae; Other 


3.81 ± 


0.44 


2.95 ± 0.45 


Clostridiurm spp. 


5.86 ± 


3.37 


0.82 ± 0.70 


Unclassified Bacteria 


2.21 ± 


0.25 


1 .84 ± 0.46 


Parabacteroides spp. 


2.25 ± 


0.12 


1 .35 ± 0.20** 


Peptostreptococcaceae; Other 


1.80 ± 


0.35 


0.93 ± 0.40 


Unclassified Firmicutes 


1.25 ± 


0.20 


0.69 ±0.13* 


Marvinbryantia spp. 


0.79 ± 


0.11 


0.81 ± 0.24 


Turicibacter spp. 


0.94 ± 


0.30 


0.52 ± 0.20 


Oscillibacter spp. 


0.62 ± 


0.29 


0.65 ± 0.28 


Asaccharobacter spp. 


0.81 ± 


0.25 


0.39 ±0.12 


Akkermansia spp. 


0.73 ± 


0.32 


0.39 ± 0.22 


Unclassified Bacteroidetes 


0.60 ± 


0.05 


0.33 ± 0.16 


Butyrici coccus spp. 


0.40 ± 


0.03 


0.52 ± 0.19 


Unclassified Alphaproteobacteria 


0.62 ± 


0.40 


0.14 ± 0.11 


Erysipelotrichaceae; Other 


0.47 ± 


0.14 


0.24 ± 0.07 


Butyricimonas spp. 


0.43 ± 


0.09 


0.22 ± 0.07 


Hyphomonadaceae; Other 


0.01 ± 


0.00 


0.58 ± 0.55 


Clostridiaceae; Other 


0.48 ± 


0.26 


0.10 ±0.08 


Alistepes spp. 


0.38 ± 


0.09 


0.18 ±0.05 


Unclassified Cyanobacteria 


0.01 ± 


0.00 


0.45 ± 0.42 


Chitinophagaceae; Other 


0.39 ± 


0.35 


0.00 ± 0.00 


Bacillariophyta; Other 


0.06 ± 


0.05 


0.28 ± 0.17 


Anaerotruncus spp. 


0.15 ± 


0.04 


0.17 ±0.02 


Coprobacillus spp. 


0.15 ± 


0.05 


0.14 ± 0.06 


Mucispirillium spp. 


0.17 ± 


0.10 


0.1 1 ± 0.05 


Ruminococcus spp. 


0.22 ± 


0.08 


0.05 ± 0.03 


Lactobacillales; Other; Other 


0.21 ± 


0.07 


0.06 ± 0.02 


Blautia spp. 


0.19 ± 


0.05 


0.07 ± 0.02 


Unclassified Bacilli 


0.18 ± 


0.04 


0.05 ± 0.02* 


Anaerostipes spp. 


0.16 ± 


0.05 


0.07 ± 0.01 


Ponticaulis spp. 


0.00 ± 


0.00 


0.22 ± 0.21** 


Anaerovorax spp. 


0.14 ± 


0.05 


0.08 ± 0.02 


Holdemania spp. 


0.16 ± 


0.06 


0.04 ± 0.01 


Bacillus spp. 


0.19 ± 


0.19 


0.01 ± 0.00 


Roseburia spp. 


0.09 ± 


0.03 


0.10 ±0.04 



Data are the mean relative abundance ± standard error. *p < .05 vs. HCC; 
^^p < .01 vs. HCC. 
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Figure 4 The absolute abundance of bacteria in thie genus Lactobacillus is reduced by exposure to the SDR stressor. Mice were exposed 
to either 0, 1, or 6 consecutive days of tine SDR stressor and qPCR used to quantify Loctobocillus group bacteria (A and B) or L reuteri (C) in tlie 
colon. *indicates p < .05 vs. 0 cycles of SDR (0 cycles is equivalent to HCC controls). Inbred C57BI/6- n = 16 at Day 0, n = 7 on Day 1, and n = 8 on 
Day 6. Outbred CD-I - n = 21 at Day 0, n = 1 5 at Day 1 , and n = 6 on Day 6. 



the majority of the fecal matter from tissue samples, 
wherein lie the majority of the luminal portion of the GI 
microbiota. This study extends previous studies that indi- 
cated that luminal and fecal microbiota can change as a 
consequence of stressor exposure to now include mucosa- 
associated populations. 

To determine whether the stressor effects only encom- 
passed changes to bacterial relative abundance or whether 
changes in absolute abundance (as assessed by determin- 
ing copies of 16 s rRNA gene per gram of sample) may 
also result from stressor exposure, qPCR was used to 
further investigate the effects of the stressor on bacter- 
ial abundance. This was performed in both an inbred 
(i.e., C57BL/6) and an outbred (i.e., CD-I) mouse strain, 
because studies consistently show that microbiota compos- 
ition is associated with mouse strain [37,38]. In addition, 
different mouse strains can also have different physiological 
and behavioral responses to stressor exposure, including 
changes to anxiety-like behavior, as well as diarrhea output 
and colonic serotonin concentration [39,40]. Thus, it 
was important to determine whether stressor effects 
were conserved across mouse strains. While 2 hrs of stres- 
sor exposure was enough to reduce the relative abundance 



Table 3 Real-time PCR assessment of bacterial group 
abundances 







C57BL6 




CD-I 




Parabacteroides distasonis 


Day 0 


6.27 (6.02 


- 7.42) 


6.65 


(6.53 - 


- 6.78) 




Day 1 


6.94 (6.66 


- 7.23) 


7.54 


(7.30 - 


- 7.77) 




Day 6 


6.52 (6.10 


- 6.94) 


7.00 


(6.26 - 


- 7.74) 


Group Porphyromonas - 


Day 0 


7.61 (7.39 


- 7.85) 


7.20 


(7.02 - 


- 7.39) 


Bocteroides - Prevotello 




8.10 (7.82 


- 8.37) 




(6.93 - 


- 7.82) 


Day 1 


7.37 




Day 6 


7.46(7.17 - 


- 7.75) 


6.97 


(6.44 - 


- 7.50) 



Data are the mean (standard error) of the copy number (loglO) of target 16 s 
rRNA per gram of tissue. 



of the genus Lactobacillus, a reduction in the absolute 
abundance was only evident after repeated exposure to 
the stressor. This indicates that some of the effects that 
stressor exposure has upon the microbiota are additive. 
Both pre- and post-stressor lactobacilli levels were similar 
in outbred CD-I mice and inbred C57BL/6 mouse strains, 
demonstrating that the effects of the stressor are consist- 
ent across mouse strains. This finding was not surprising 
given that stressor-induced reduction in lactobacilli levels 
have been documented in other host species, including 
human and non-human primates. 

The finding that stressor exposure reduced the relative 
and absolute abundance of lactobacilli is consistent with 
previous studies demonstrating that stressor exposure 
can reduce the lactobacilli [19,31,35,36]. To date, however, 
the Lactobacillus species that can be reduced by stressor 
exposure has not been addressed. Many species of bacteria 
in the genus Lactobacillus are known to have health pro- 
moting effects, and studies from this laboratory, as well as 
others, indicate that administering probiotic L, reuteri to 
mice can reduce colonic inflammatory responses [20,21,41]. 
Because stressor exposure can exacerbate colonic inflam- 
mation [42], we focused on whether indigenous L, reuteri 
was reduced in stressor-exposed animals. As predicted, L, 
reuteri was significantly reduced in colonic mucosa from 
outbred CD-I mice exposed to the SDR stressor. Interest- 
ingly, however, L, reuteri was not consistently detectable 
in inbred C57BL/6 mice. Thus, it is apparent that L, reu- 
teri levels on the mucosa of C57BL/6 mice are lower than 
levels observed in outbred CD-I mice and below the de- 
tection limit (4.5 log copies/gram of wet tissue). 

Although it is difficult to compare deep sequencing re- 
sults across different experiments that assess different 
body niches and encompass different stressors, the use 
of 454 pyrosequencing demonstrated that the effects of 
stressor exposure are not confined to effects on the 
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Table 4 Real-time PCR assessment of colonic inflammation 



C57BL6 CD-I 







DAY 1 




DAY 6 


DAY 1 


DAY 6 


ILip 


HCC 


1.00 (0.57 - 


1.75) 


1.00 (0.55 - 1.82) 


1 .00 (0.37 - 2.68) 


1.00 (0.64 - 1.57) 




SDR 


0.95 (0.49 - 


1.84) 


1.37 (1.00 - 1.37) 


3.18 (2.92 - 3.47) 


1.11 (0.59 - 2.11) 


iNOS 


HCC 


1 .00 (0.56 - 


1.78) 


1.00 (0.50 - 2.01) 


1 .00 (0.56 - 1 .78) 


1 .00 (0.69 - 1 .46) 




SDR 


0.62 (0.46 - 


0.83) 


2.50 (1.92 - 3.25) 


0.32 (0.24 - 0.43) 


0.23 (0.14 - 0.37) 


TNFa 


HCC 


1 .00 (0.74 - 


1.35) 


1.00 (0.51 - 1.98) 


1 .00 (0.76 - 1 .32) 


1.00 (0.61 - 1.65) 




SDR 


0.50 (0.35 - 


0.73) 


1.03 (0.71 - 1.49) 


0.92 (0.69 - 1 .22) 


0.92 (0.63 - 1 .33) 



Data are the mean (standard error) of the fold change in cytokine gene expression over non-stressed HCC controls. 



lactobacilli. The relative abundance of bacteria in the 
family Porphyromonadaceae was found to be significantly 
reduced by exposure to the SDR stressor. This finding is 
consistent with a previous study demonstrating that ex- 
posure to a prolonged restraint stressor in oubred CD-I 
mice was sufficient to reduce the relative abundance of 
cecal Porphyromonadaceae [23]. However, it was not 
previously known whether stressor exposure could re- 
duce the absolute abundance of the Porphyromonada- 
ceae family. Thus, qPCR was used to determine whether 
this stressor also affected the absolute abundance of 
Porphyromonadaceae by assessing levels of bacteria in 
the Porphyromonas-BacteroideS'Prevotella group as a 
surrogate marker. This was performed because of the 
problems designing qPCR primers specific for the Por- 
phyromonadaceae family. The absolute abundance of the 
Porphyromonas-BacteorideS'Prevotella group was similar 
in CD-I mice and inbred C57BL/6 mice, and stressor ex- 
posure had no effect on the absolute abundance in either 
mouse strain. Thus, although this study and others [23] 
have found stressor-induced reductions in the relative 
abundance of Porphyromonadaceae, these results indicate 
that absolute abundance of the Porphyromonadaceae is 
unaffected by stressor exposure. Stressor exposure also 
did not affect the absolute abundance of P. diastonosis 
in either strain of mice, even though the relative abun- 
dance of Parabacteroides was reduced as detected by 
454 pyrosequencing. 

Colonic inflammation is well known to impact the gut 
microbiota, and studies have found that intestinal inflam- 
mation can also reduce the abundance of the family Lacto- 
bacillaceae [43,44] and Porphyromonadaceae [45]. Given 
the bidirectional interaction between the microbiota and 
intestinal inflammation and increased intestinal physio- 
logical inflammation upon exposure to stressor [46], it is 
often suggested that increases in intestinal inflammation 
during stressful periods could be responsible for changes 
in gut microbiota. However, inflammatory cytokine gene 
expression was not significantly changed by either 2 hrs 
of the SDR stressor or 2 hrs of the stressor repeated on 
6 consecutive days. Thus, the results of this study 



indicate that it is not likely that stressor-induced colonic 
inflammation is responsible for altering the abundances 
of the indigenous colonic mucosa-associated microbiota. 

The effects of SDR on male mouse immunological, 
physiological, and behavioral functioning have been well 
characterized over the past fifteen years [47-49]. The con- 
sistent activation of the HPA axis and the SNS makes SDR 
an ideal method of inducing an acute stress response in 
male mice. While it would be desirable to determine the 
effects on female mice as well, aggression in female mice is 
too low to induce a physiological stress response in this 
paradigm. Other stressor paradigms, however, have shown 
that stressor exposure in a lab setting can also affect the 
microbiota composition of female mice [50,51]. Thus, it is 
likely that a social stressor in female mice would also 
impact the lactobacilli, as well as Porphyromonadaceae and 
Parabacteroides, Interestingly, studies have shown that 
male and female microbiota can evoke different levels of 
host sex hormone release and immunological outputs that 
can then feedback and alter microbial composition [52]. 
Future studies comparing the microbiota, mucosal immu- 
nity and endocrine responses in both male and female mice 
are likely to demonstrate additional multidirectional inter- 
actions between host physiology and the microbiota. This 
requires further study. 

Conclusions 

Gut microbes have important effects on the host espe- 
cially when they are associated with the mucosa. Previous 
studies have found that stressor exposure can affect the 
microbiota, in particular the lactobacilli. However, many 
of these studies focused on microbes that are shed in the 
stool during prolonged or chronic stressors [19,23,31,50]. 
It was not previously clear whether the effects of stress on 
the microbiota were limited to changes in the number of 
bacteria shed in the stool, nor whether short lasting acute 
stressors could also impact the microbiota. The present 
study demonstrates that exposure to as little as 2 hrs of a 
social stressor is sufficient to significantly affect some pop- 
ulations of the colonic mucosa-associated microbiota. 
This builds upon previous work by showing that exposure 



Galley et al. BMC Microbiology 2014, 14:189 
http://www.bionnedcentral.conn/1471 -21 80/1 4/1 89 



Page 8 of 1 3 



to the social stressor affects multiple regions (e.g. colon, 
cecum) and niches (e.g. lumen, mucosa-associated) of the 
GI tract, as well as two different strains of mice, suggest- 
ing that the effect of social stresss exposure upon the 
microbiota is a universal process, as opposed to either an 
artifact or isolated finding. The mechanisms by which this 
occurs are not yet understood, but are not likely to involve 
stressor-induced inflammation in the colon. Future studies 
should assess the impact of stressor- associated hormones 
and gut functioning since others have demonstrated that 
activation of the stress-associated hormone corticotrophin 
releasing hormone (CRH) alters gut microbiota, an effect 
that was associated with alterations in intestinal motility 
[53]. Because the gut microbiota has been associated with 
diverse diseases that are known to be exacerbated by 
stressful situations, such as the inflammatory bowel dis- 
eases, irritable bowel syndrome, and even multiple scler- 
osis, it is becoming increasingly important to understand 
the impact that stress can have on the microbiota and 
whether stressor-induced alterations in the microbiota are 
involved in stressor-induced disease exacerbation [54-57]. 

Methods 

Mice 

Male CD-I mice and male C57BL/6, aged 6-8 weeks, were 
ordered from Charles River Laboratories (Wilmington, 
MA; Raleigh, NC). The mice were housed in groups of 
2 or 3 mice per cage in an approved vivarium and 
were allowed to habituate to the vivarium for one 
week prior to testing. The cages were kept in an approved 
vivarium with a 12:12 hour light-dark schedule with lights 
on 0600 to 1800. Food and water was available ad libitum. 
The Ohio State University's Animal Care and Use Com- 
mittee approved all experimental procedures (Protocol 
2009A0235-R1). CD-I mice were used in cytokine and 
bacterial qPCR experiments. C57BL/6 mice were only 
used in the 454 pyrosequencing analysis and the cytokine 
and bacterial qPCR experiments. 

Social disruption stressor 

The SDR experiments were performed as described pre- 
viously [58]. Briefly, an aggressive male CD-I retired 
breeder, termed the aggressor, was placed into a cage 
with younger C57BL/6 resident mice at 1700 hours, 
which represents the beginning of the mouse active 
cycle. The aggressor repeatedly attacked and defeated 
the C57BL/6 test mice over the course of the 2 hour 
SDR cycle. If the aggressor did not begin to attack the 
test mice within 10 min of placement in the cage, it was 
removed and another aggressor was put in. Wounding 
was monitored throughout the SDR cycle. Only slight 
superficial wounds were allowed. If wounds that pene- 
trated the cutaneous layer developed over the course of 
SDR, those test mice were removed from the study. While 



it would be desirable to test the effects of the stressor 
on both male and female mice, female mice do not fight 
in this paradigm. Thus, all analyses were limited to male 
mice. 

Tissue removal 

After exposure to the SDR stressor (a single 2 hr cycle 
for pyrosequencing, a single 2 hr cycle or 2 hr cycles re- 
peated on 6 consecutive days for qPCR), or at the equiva- 
lent time in non-stressed HCC mice, mice were euthanized 
using CO2 asphyxiation. Colons were removed aseptically. 
Luminal contents were removed by cutting into the tissue 
where liquid and/or solid contents rested, and gently 
extracting with forceps. This represented the non-mucosa 
associated populations. Scraping/tissue squeezing was not 
performed so as not to disturb the mucosa-associated pop- 
ulations, which are adhered to an easily disturbed mucous 
layer. The colonic tissue and remaining adherent micro- 
biota populations were placed into a microcentrifiige tube 
and snap frozen in liquid nitrogen. All samples were stored 
at -80°C until processing. Studies involving 454 pyrotag se- 
quencing utilized the C57BL/6 colonic tissue that was col- 
lected after a single 2 hr cycle of SDR, with n = 5 per group 
(HCC vs. SDR). Animals were kept 3 and 2 mice/cage to 
minimize cage effects. Studies involving bacterial quantifi- 
cation PCR utilized colonic tissue collected at 0 days of 
SDR (n=16 C57BL6, n = 21 CD-I), following a single 
2 hr cycle of SDR (n = 7 C57BL/6, n = 15 CD-I) as well as 
samples collected after SDR repeated on 6 consecutive days 
(n = 8 C57BL/6, n = 6 CD-I). All experiments were n = 3 or 
4 per group. The single 2 hour cycle of SDR in C57BL/6 
was replicated once, and six repeated exposures of SDR 
on C57BL/6 was replicated twice. The single 2 hour cycle 
of SDR in CDl mice was replicated four times, and six re- 
peated exposures of SDR on CDl mice was replicated 
once. Successive replicates were performed to reduce high 
variability between samples in bacterial qPCR. No-SDR- 
exposure groups were sacrificed alongside both single 
2 hr cycle of SDR and 6 repeated cycles of SDR groups 
and combined during statistical analyses, resulting in an 
increased number of replications over test groups. 

bTEFAP pyrosequencing 

Amplicon pyrosequencing (bTEFAP®) was originally 
described by Dowd et al. on C57BL/6 mice colonic 
samples and has been utilized in describing a wide range 
of environmental and health related microbiomes includ- 
ing the intestinal populations of a variety of sample types 
and environments [59-61]. In this protocol, a 1-mm 
segment was used from the center of the entire colonic 
length. In a modified version of this process, 16S universal 
Eubacterial primers 530 F 5'GTGCCAGCMGCNGCGG 
and llOOR 5'GGGTTNCGNTCGTTR were used in a 
single-step 30 cycle PCR using HotStarTaq Plus Master 
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Mix Kit (Qiagen, Valencia, CA). The following conditions 
were used: 94°C for 3 minutes, followed by 28 cycles of 
94°C for 30 seconds; 53°C for 40 seconds and 72°C for 

1 minute; after which a final elongation step at 72°C for 
5 minutes was performed. Following PGR, all amplicon 
products from different samples were mixed in equal con- 
centrations and purified using Agencourt Ampure beads 
(Agencourt Bioscience Corporation, MA, USA). Samples 
were sequenced utilizing Roche 454 FLX titanium instru- 
ments and reagents and following manufacturer s guidelines. 

Microbial community analysis 

The returned sequences were processed and analyzed using 
Quantitative Insights Into Microbial Ecology (QIIME) 
(version 1.4.0) [62] according to standard protocols. Filter- 
ing passed sequences based on: length between 200 and 
1000 bp, zero ambiguous bases or primer mismatches, 
and no homopolymer runs greater than 6. 97.9% of the se- 
quences passed quality filtering for further analysis. After 
filtering, there was an average of 10496 sequences per 
sample. Sequences with 97% similarity were grouped into 
OTU using UCLUST [63]. Representative sequences from 
each OTU were aligned against the Greengenes core refer- 
ence alignment using PyNAST [64,65]. The RDP Classifier 
was used to assign taxonomy to each representative se- 
quence against the RDP database using standard options 
[66]. The minimum confidence threshold for taxonomic 
assignment was 0.80. OTUs were considered unclassified 
if there was not a strong match within this confidence 
interval for the representative sequence within the RDP 
database. A phylogenetic tree was built using FastTree 
from the aligned OTU representative sequences for deter- 
mining UniFrac distances between samples [67]. 

Quantification of bacterial groups using quantitative 
polymerase chain reaction (qPCR) 
DNA extraction 

50 mg of tissue, taken from the ascending colon portion 
of the colonic tract, were homogenized in liquid nitrogen 
in a mortar and pestle and transferred to a microcentri- 
fuge tube. The samples were incubated for 75 min in lyso- 
zyme buffer (20 mg/mL lysozyme, 20 mM Tris-HCl, 

2 mM EDTA, 1.2% Triton) at 37°C to weaken the cell wall 
of Gram-positive bacteria. Following this, DNA was iso- 
lated using a modified version of the Qiagen DNA Mini 
Isolation Kit protocol (Qiagen, Germany). In short, 100 [iL 
of Buffer ATL and 20 [iL of proteinase K were added to 
the microcentrifuge tube and incubated for 3.5 hr at 56°C, 
with manual vortexing at 30 minute intervals for the dur- 
ation of the incubation. The samples were heated at 85°C 
for 10 min to inactivate proteinase K and then transferred 
into a new tube containing 20 [iL of proteinase K and 
200 [iL of Buffer AL. The samples were incubated at 56°C 
for 30 min followed by 95°C for 10 min. 200 [iL of 



absolute ethanol was then added to the sample and the 
sample was added to a QIAamp spin column. After 
sample washing, the DNA from each sample was recov- 
ered into 100 [iL of Buffer AE using a QIAamp spin col- 
umn per manufacturers protocols. The DNA extracts 
were stored at -80°C until analysis was performed. Gen- 
omic DNA from pure bacterial cultures was isolated es- 
sentially as described above, but using a Qiagen 
protocol for the isolation of genomic DNA from Gram- 
positive bacteria and shorter (45 minutes instead of 
75 minutes) incubation in the lysozyme buffer. 

Sample derived standards 

The method of Chen et al. [68] was adapted for the gener- 
ation of sample-derived qPCR standards [68]. Briefly, 
equivalent mass amounts of colonic DNA from each sam- 
ple were pooled as the template for PGR to obtain even 
bacterial representation in order to generate standards for 
bacterial copy number quantification. DNA was quantified 
using the Quant-it PicoGreen kit (Life Technologies, 
Grand Island, New York). Amplification PGR targeting 
specific genera or species-level 16S rRNA was performed 
to obtain standard amplicons based upon the representa- 
tive populations within each sample. PGR was achieved 
using an ABIPrism 7000 thermocycler. The primer sets 
are shown in Table 5. The standard amplification PGR 
thermoprofile for the genus Lactobacillus consisted of ini- 
tial denaturation at 94°G for 5 min, followed by 35 cycles 
of 94°G for 30 sec, 58°G for 30 sec, and 72°G for 40 sec, 
and final extension at 72°G for 10 min. For L. reuteri PGR, 
the thermoprofile used was the same as those for the 
genus-specific except for a higher primer annealing 



Table 5 PCR primers and probes 



Lactobacillus genus 


Forward 


AGCAGTAGGGAATOTCCA 




Reverse 


CACCGCTACACATGGAG 


Lactobacillus reuteri 


Forward 


CAGACAATCmGATOmAG 




Reverse 


GOTGTOGmGGGCTOTC 


Parabacteroides distasonis 


Forward 


TGCCTATCAGAGGGGGATAAC 




Reverse 


GCAAATATOCCATGCGGGAT 


Porphyromonas-Bacteroides- 


Forward 


GGTGTCGGCTOAGTGCCAT 


Prevotella 


Reverse 


CGGA(C^GTAAGGGCCGTGC 


TNF-a 


Forward 


CTGTCTACTGAAOTCGGGGTGAT 




Reverse 


GCTCTGGGCCATAGAACTGATG 




Probe 


ATGAGAAGTTCCCAAATGGCCTCCCTC 


IL-ip 


Forward 


GGCCTCAAAGGAAAGAATCTATACC 




Reverse 


GTATOOTGGGATCCACACTCT 




Probe 


ATGAAAGACGGCACACCCACCCTG 


iNOS 


Forward 


CAGCTGGGCTGTACAAACOT 




Reverse 


TGAATGTGATGmGCTOGG 




Probe 


CGGGCAGCCTGTGAGACCmGA 
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temperature (60°C instead of 58°C). The Bacteroides- 
Prevotella-Porphyromonas standard amplification thermo- 
profile was 1 cycle at 95°C for 5 min, followed by 35 cycles 
of 95°C for 20 s, 68°C for 30 s, and 72° for 55 s, with 1 
final annealing cycle of 72°C for 5 min. The Parabacter- 
oides distasonis standard amplification thermoprofile was 

1 cycle at 94°C for 5 min, followed by 30 cycles of 94°C for 
30 s, 60°C for 30 s, and 72°C for 40 s, with 1 final anneal- 
ing cycle at 72°C for 8 min. After the amplification PGR, 
the PGR product was purified using the Qiagen PGR 
Purification Kit and quantified by Quant-it PicoGreen kit 
to determine copy number of the community standards. 

A standard curve for absolute quantification was cre- 
ated from the previously amplified standards to encom- 
pass 10^-10^ copies per reaction. The population of total 
lactobacilli and L, reuteri was separately quantified using 
respective specific primers and SYBR Green against 
these dilution series as described previously [69]. qPGR 
thermoprofiles for Lactobacillus and L, reuteri were 1 cycle 
at 94°G for 5 min, followed by 40 cycles of 94°G for 30 s, 
specific annealing temperature for 30 s, and 72°G for 40 s. 
Annealing temperature for Lactobacillus was 58°G and L, 
reuteri was 60°G. qPGR thermoprofile for P, distasonis 
was 1 cycle at 50°G for 2 min, 1 cycle at 95°G for 10 min, 
then 40 cycles at 95°G for 15 s and 60°G for 60 s. qPGR 
thermoprofile for Bacteroides-Prevotella-Porphyromonas 
was 1 cycle at 95°G for 5 min, followed by 40 cycles at 95°G 
for 15 s, 68°G for 20 s, and 72°G for 30 s. Abundance of 
bacterial groups (copies of 16 s rRNA gene per gram/ 
sample) was computed based on the copies of qPGR reac- 
tion and the number of reactions that could be performed 
with the DNA derived from 1 g of each tissue sample. 
The detection limit were as follows: total lactobacilli 
and Bacteroides-Prevotella-Porphyromonas was -4.0 log 
copies/gram of wet tissue (logio), L, reuteri was -4.5 log 
copies/gram of wet tissue (logio), P- distasonis was -5.0 
log copies/gram of wet tissue. Primer sequences have been 
previously published and are listed in Table 5 [70-74]. 

Quantification of colonic cytokine and inflammatory 
mediator mRNA using quantitative RT-PCR 

Total RNA was isolated from whole colonic tissue samples 
(-120 mg) using TRI-zol (Invitrogen, Garlsbad, Galifornia) 
according to manufacturers instructions. The extracted 
RNA was quantified via spectrophotometry, and con- 
verted to cDNA using the Reverse Transcription System 
(Promega, Madison, WI). qPGR was then completed using 
a master mix containing 2x Universal TaqMan master 
mix (Life Technologies, Grand Island, New York), 0.9 (iM 
(each) forward and reverse primers (see Table 5), and 
0.250 (ig sample cDNA. 18S was used as the housekeeping 
gene. The PGR was performed using a Prism 7000 se- 
quence detection system with the following thermoprofile: 

2 min at 50°G, 10 min at 95°G, and then 40 amplification 



cycles of 15 s at 95°G and 1 min at 60°G. The relative 
amount of mRNA was determined using the comparative 
cycle threshold (C^) method as previously described [75]. 
Non-stressed HGG control samples were used as baseline 
controls and were set at a value of 1. All other samples are 
based on a fold change from these control samples. Pri- 
mer sequences have been previously published and are 
listed in Table 5 [76] . All groups were n = 3 - 5, using a 
single experiment without replicates. 

Statistical analyses 

Alpha diversity measurements, including Shannon diver- 
sity index, equitability (evenness) and Ghao (richness) 
were computed using QIIME on the groups (HGG control 
vs. SDR Stressor) [77]. Non-parametric T- tests with 999 
Monte Garlo permutations were used on these alpha di- 
versity measurements to test for significance at a sampling 
depth of 4805 sequences per sample. 

Taxonomic abundances attained from QIIME were 
compared between the two groups of mice using non- 
parametric Mann- Whitney U Tests. PGoA of unweighted 
UniFrac distance matrices was used to determine clustering 
between the two groups (HGG control vs. SDR Stressor) 
[78]. Analysis of Similarity (ANOSIM), a beta-diversity stat- 
istic that is available through the vegan package of R and 
accessible with QIIME, was used to calculate statistical 
significance between the distance matrices of groups at 
999 permutations [79,80]. All of these analyses were per- 
formed in QIIME. 

Differences in bacterial abundances were determined 
using non-parametric Kruskal Wallis tests with the cycle 
of SDR (i.e., 0, 1, or 6 cycles) as the single factor. Mann- 
Whitney U tests were performed a priori, Ghanges in 
cytokine gene expression were determined using a two- 
factor ANOVA with group (HGG control vs. SDR Stressor) 
and Day (1 Day vs. 6 Days) as the two independent factors. 
Protected means comparisons were used as post-hoc tests. 
In all cases, a was set at 0.05, while tendency was declared 
at 0.05 < p < 0.10. These tests were performed using SPSS 
V.21 (IBM, Ghicago, IL). 

Availability of supporting data 

The sequences supporting the results of this article are 
available in the NGBI Sequence Read Archive under the 
study accession number SRP035598 (http://www.ncbi. 
nlm.nih.gov/sra/?term=SRP035598). 

Abbrevations 

Gl: Gastrointestinal; SDR: Social disruption; qPCR: Quantitative real-time 
polymerase chain reaction; HPA: Hypothalamic-pituitary-adrenal axis; 
OTU: Operational taxonomic unit; UPGMA: Unweighted pair-group method 
with arithmetic mean; ANOSIM: Analysis of similarity; HCC: Home cage 
control; IL-lp: Interleukin 1-beta; TNF-a: Tumor necrosis factor-alpha; 
iNOS: Inducible nitric oxide synthase; DSS: Dextran sodium sulfate; 
CRH: Corticotrophin-releasing hormone; bTEFAP: Bacterial tag-encoded FIX 
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amplicon pyrosequencing; QIIME: Quantitative Insights Into Microbial 
Ecology; PCoA: Principle coordinates analysis. 
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